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ABSTRACT 

The detection of narrow SiO thermal emission toward young outflows has been proposed to be a 
signature of the magnetic precursor of C-shocks. Recent modeling of the SiO emission across C- 
shocks predicts variations in the SiO line intensity and line shape at the precursor and intermediate- 
velocity regimes in only few years. We present high-angular resolution (3.8" x3.3") images of the 
thermal SiO J=2— )-l emission toward the L1448-mm outflow in two epochs (November 2004-February 
2005, March-April 2009). Several SiO condensations have appeared at intermediate velocities (20- 
50kms~^) toward the red-shifted lobe of the outflow since 2005. Toward one of the condensations 
(clump D), systematic differences of the dirty beams between 2005 and 2009 could be responsible for 
the SiO variability. At higher velocities (50-80 kms~^), SiO could also have experienced changes in 
its intensity. We propose that the SiO variability toward L1448-mm is due to a real SiO enhancement 
by young C-shocks at the internal working surface between the jet and the ambient gas. For the 
precursor regime (5.2-9.2 kms~^), several narrow and faint SiO components are detected. Narrow SiO 
tends to be compact, transient and shows elongated (bow-shock) morphologies perpendicular to the 
jet. We speculate that these features are associated with the precursor of C-shocks appearing at the 
interface of the new SiO components seen at intermediate velocities. 

Subject headings: star: formation — physical data and processes: shock waves — ISM: jets and 
outflows — ISM: individual (L1448) 



1. INTRODUCTION 

Variability is a common phenomena in the evolution 
of very young stars. Joy (1945) very early found that 
T Tauri stars have an irregular behavior in their op- 
tical emission. More recently, simultaneous X-ray and 
optical/near-IR variability studies have shown that a 
high fraction of pre-main-sequence stars are variable in 
X-ray and optical/near-IR emission (Stassun et al. 2006; 
Forbrich et al. 2007). At radio wavelengths, young stel- 
lar objects are known to show variability in their radio 
continuum emission in time-scales of few years (Curicl 
et al. 2006; Pech et al. 2010). In the case of molecular 
line emission, only non-thermal maser emission has been 
found to be variable (see Reid & Moran 1981; Elitzur 
1992). However, variability in the morphology, line in- 
tensity and line profiles of thermal molecular emission 
still remains to be reported. 

Broad SiO thermal line emission is known to be a good 
indicator of shocked gas in molecular outflows (Martin- 
Pintado et al. 1992). The SiO line profiles have typi- 
cal linewidths of tens of kms~^ and appear shifted with 
respect to the ambient cloud velocity. In addition to 
this broad emission, Lefloch et al. (1998) and Jimenez- 
Serra et al. (2004) reported the detection of very narrow 
SiO thermal lines (linewidths of <lkms~"'^) centered at 
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almost ambient velocities toward the young NGC1333 
and L1448-mm outflows. Toward L1448-mm, Jimenez- 
Serra et al. (2004) proposed that narrow SiO could arise 
from the magnetic precursor region within very young 
C-shocks (Jimenez-Serra et al. 2004). 

At the early stages of the C-shock propagation, the 
magnetic precursor of C-shocks accelerates, compresses 
and heats the ion fluid before the neutral one (Draine 
1980). The subsequent ion- neutral velocity drift gener- 
ates enough friction to sputter dust grains (Caselli et al. 
1997; Gusdorf et al. 2008), progressively enhancing the 
gas-phase abundances of SiO (Jimenez-Serra et al. 2008). 
Jimenez-Serra et al. (2009) have claimed that the SiO line 
proflles in young outflows should then evolve from very 
narrow lines arising from the magnetic precursor region, 
to broader SiO emission generated in the intermediate- 
and high- velocity postshock gas. From this, variability of 
the narrow SiO lines and of the intermediate- velocity SiO 
emission should appear in only few years (Jimenez-Serra 
et al. 2009). 

We present interferometric images of the SiO J=2^>1 
thermal line emission toward the L1448-mm outflow ob- 
served with the PdBI at two different epochs (2005 and 
2009)^. These images show that several SiO condensa- 
tions (named A, B, C, D and E) have recently appeared 
at intermediate velocities along the direction of the red- 
shifted lobe of the outflow since 2005. The increase in the 
SiO flux toward condensations A, B and C are likely real 
and produced by an enhancement in the gas-phase abun- 
dance of SiO by the recent interaction of young C-shocks 
with the ambient gas. The variable nature of the SiO 
clumps D and E, however, remains unclear (Section 4.1). 

'^ Based on observations carried out with the IRAM Plateau 
de Bure Interferometer. IRAM is supported by INSU/CNRS 
(France), MPG (Germany) and IGN (Spain) 
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Our results also reveal the detection in 2009 of several 
new narrow SiO components that tend to be found at 
the interface between the SiO condensations A, B, C, 
D and E reported at intermediate velocities toward the 
red-shifted lobe of the outflow (Section4.2). We specu- 
late that the recently detected narrow SiO components 
are possibly associated with the magnetic precursor of C- 
shocks (Sections). This is the first time that variability 
of molecular thermal emission is ever reported. 

2. OBSERVATIONS 

PdBI observations of the thermal SiO v—0 J=2— >1 
line emission (at 86.846 GHz) toward the L1448-mm out- 
flow, were carried out in the B and C configurations 
in November 2004-February 2005 (for ~15hrs of on- 
source time) and in the C configuration in March- April 
2009 (^9hrs). The phase center of the images was set 
at a(J2000)=03''25'"38.8^ (5(J2000)=-h30°44'05". The 
correlator setup included two different spectral resolution 
units (with bandwidths of 20 MHz and 80 MHz), which 
provided spectral resolutions of 39 kHz and 312 kHz, re- 
spectively (i.e. velocity resolutions of 0.14 and 1.1 kms~^ 
at ^87 GHz). We used 3C84 (9.5 Jy) as bandpass calibra- 
tor, and MWC349 (1.2 Jy) as flux calibrator. 0333+321 
(2 Jy) and 0234-f-285 (2 Jy) were used as gain calibra- 
tors. Pointing and focus of the antennas was checked 
every 65 min and re-adjusted if necessary. Pointing er- 
rors were typically ^2" (i.e. less than 4% the field-of- 
view, FOV, of 58" at -87 GHz) in the 2005 observations, 
and ~2-4" (-4-8% the FOV) in the 2009 tracks. Focus 
corrections were <0.1mm in the 2005 observations, and 
—0.05-0.15 mm in the 2009 tracks. Data calibration was 
done with the GILDAS package (CLIC)"^. We used the 
"Data Quality Assessment" option available within this 
software, to make sure that only data with a phase noise 
smaller than 40° were retained for constructing the final 
maps. 

3. UV DATA ANALYSIS AND IMAGING 

Construction of the UV tables, data imaging, con- 
tinuum subtraction and cleaning were also carried out 
within GILDAS by using CLIC and MAPPING. Inde- 
pendent UV tables were generated for the 20 MHz and 
80 MHz correlator units used in the observations. The 
UV data were resampled in velocity by imposing the 
same number of channels (460 and 230 for the 20 MHz 
and 80 MHz units, respectively), the same position for 
the central channel (230 for the 20 MHz unit, and 115 for 
the 80 MHz unit), LSR velocity of the source (— 5kms~^ 
in all cases) and channel spacing (39kHz or 0.14 km s~^ 
for the 20 MHz unit, and 312kHz or 1.1 km s^^ for the 
80 MHz unit). 

3.1. The 80 MHz low spectral resolution data 

The SiO emission toward the L1448-mm outfiow is 
known to be detected to velocities up to ±80kms~^ 
(Martin-Pintado et al. 1992; Jimenez-Serra et al. 2004, 
2005). The velocity coverage provided by the 80 MHz 
correlator units was — 276kms~^, covering the entire 
SiO line profile at low spectral resolution (312 kHz or 
l.lkms"^). The morphology of the SiO J=2^-l emis- 
sion at intermediate (from 20 to 50kms^^) and high ve- 
locities (from 50 to 80kms^^) is extended (Guilloteau et 



al. 1992; Girart & Acord 2001) and therefore sensitive to 
variations of the UV plane coverage. In particular, since 
the configuration of the PdBI array in 2009 was more 
compact than in 2005, the 2009 observations could have 
been more sensitive to extended emission introducing a 
low-level contribution that could have mimicked a time- 
variation in the flux and morphology of the SiO emission. 
To avoid this potential problem, we have homogenized 
the 2005 and 2009 datasets in the UV plane by selecting 
the common visibilities to both epochs contained within 
UV ceUs of 10 m-size. We created a MAPPING proce- 
dure that sets the weight of the non-overlapping visi- 
bilities to negative values so that these visibilities were 
flagged out and not considered in the subsequent steps 
for imaging. As shown below, this has the disadvantage 
to decrease the sensitivity of the original observations, 
but guarantees that any possible variation detected in 
the SiO J=2— )-l line profile toward L1448-mm is not af- 
fected by differential missing flux due to different UV 
sampling in the two epochs. 

The continuum UV tables were generated by consider- 
ing the same line-free frequency ranges in both datasets. 
After selecting the common visibilities in both contin- 
uum UV tables, we subtracted the continuum UV data 
from the line+continuum UV tables for the 2005 and 
2009 epochs by using the UVSUBTRACT task within 
MAPPING. 

The dirty images were created using the UVMAP 
task within MAPPING. The size of the dirty maps 
was 256 pixelsx 256 pixels, with pixel size of 0.86" xO.86". 
We used robust weighting to supress the contribu- 
tion from the secondary lobes of the beams. The re- 
sulting synthesized beams were very similar (see also 
Section4.1). For the line dirty images, the beam 
sizes were of 3.80"x3.35", P.A.=179° for the 2005 im- 
ages, and 3.87"x3.31", P.A.=157° for the 2009 maps. 
For the continuum datacubes, the synthesized beams 
were 3.48" x 3.33", P.A.= 28° for the 2005 epoch, and 
3.44" X 3.37", P.A.=129° for the 2009 images. 

We finally cleaned the dirty images by using the Hog- 
bom method within MAPPING. We considered for all 
images a circular support of 60" -diameter (comparable 
to the FOV) . The stopping criteria for the CLEAN task 
was met for every channel when the maximum amplitude 
of the absolute value of the residual image was lower than 
2.5% the peak intensity of the dirty map. We note that 
this criterion is set up by default within the MAPPING 
software. The final analysis of the cleaned images was 
carried out with the new MADCUBAJ software*. 

An automatic data analysis procedure was applied to 
the UV data to calculate the visibility noise for every 
single channel in the UV tables. In Table 1, we show 
the median values obtained from the histograms of the 
visibility noise per channel for the 80 MHz unit tables. 
The noise is almost a factor of 2 larger in the 2005 
UV data (— 23mJy/beam per 312 kHz- channel) with re- 
spect to the 2009 dataset (— 13mJy/beam per 312 kHz- 
channel). Hereafter, we will use the thermal noise of the 
2005 data as a reference to search for differences in the 
morphology and intensity of the SiO line emission at in- 
termediate and high velocities between both epochs. We 
note that the thermal noise of the original 2005 and 2009 



'' See http://www.iram.fr/IRAMFR/GILDAS. 
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UV tables (i.e. considering the common and the non- 
overlapping UV data) ranged between 5 and 7 mJy/beam 
per 312kHz-channel, i.e. factors of 2-3 better than those 
obtained in the common UV data tables. 

Masson (1986) proposed another technique to mea- 
sure very small changes in the structure of very bright 
radiocontinuum sources such as the planetary nebula 
NGC7027. The cross-calibration method is based on the 
use of the closure phases and amplitudes derived from 
the data from epoch I, to minimize the effects of calibra- 
tion errors on the final subtracted map. This method, 
however, assumes that the source is bright enough to do 
self-calibration. The 3 mm continuum flux measured to- 
ward the L1448-mm source is only ~10mJy, insufficient 
to apply this technique (see Section4.1). 

3.2. The 20 MHz high spectral resolution data 

With the high-spectral resolution 20 MHz units, we im- 
aged the very narrow SiO J=2— >1 line emission found 
toward the L1448-mm outfiow at velocities close to the 
ambient cloud velocity (Jimenez-Serra et al. 2004). This 
narrow emission is expected to be faint and therefore, the 
data imaging procedure followed in Section 3.1 cannot be 
applied to the 20 MHz data because of the poor sensitiv- 
ity achieved in the final images. We believe that differ- 
ences in the UV coverage between both epochs do not 
represent any problem for our variability studies of the 
narrow SiO lines toward L1448-mm, because this emis- 
sion is found to be compact and to arise from practically 
unresolved condensations associated with the red-shifted 
lobe of the outflow (see Section4.2). In any case, we only 
used the C configuration 2005 UV data so that the UV 
coverage in the 2005 images was kept as close as possible 
to that obtained in the 2009 observations (also carried 
out in C configuration). 

The 2005 and 2009 datasets were subsequently ho- 
mogenized by weighting the 2005 UV-data (beam of 
3.52"x2.22",' P.A.=35°) with a Gaussian taper of 
80 mx 100m and at a P.A.=65°. The resulting 2005- 
beam was 5.20" x 3.65", P.A.=21°, which matched the 
2009-beam of 5.18"x3.69", P.A.==57°. The map size 
and pixel size used to construct the dirty images were 
the same as those of the 80 MHz data images. In this 
case, however, we used natural weighting instead of ro- 
bust weighting to increase the final signal-to-noise ratio. 
The dirty images were cleaned in a similar fashion as the 
80 MHz maps. The derived visibility noise per channel 
for the 20 MHz UV datasets is shown in Table 1. Like in 
the 80 MHz data, the noise in the 2005 observations was 
higher (21.8 mJy/beam per 39 kHz-channel) than in 2009 
(18.7mJy/beam per 39kHz-channcl). We will then use 
the 2005 thermal noise level as the reference to search for 
variability in the narrow SiO lines toward L1448-mm. 

We note that, although the noise in the high-spectral 
resolution 20 MHz data is expected to be significantly 
worse than in the 80 MHz UV data, the 20 MHz data 
considers all visibilities available from the 2005 and 2009 
observations, while the low-spectral resolution 80 MHz 
UV files have been constructed from only UV data com- 
mon to both epochs. 

4. RESULTS 

4.1. Variability of the thermal SiO emission at the 
intermediate- and high-velocity shock regimes 



In Figure 1 , we present the averaged intensity images 
of the SiO J=2^1 emission from 20 to 50kms^^ ob- 
served in 2005 (left panel) and 2009 (central panel) 
toward the red-shifted lobe of the L1448-mm outflow. 
For our analysis, we use averaged intensity images in- 
stead of integrated intensity maps because it will al- 
low a simpler comparison of the data with the response 
of the interferometer (beam) for both epochs (see be- 
low). Contour levels are the same in both averaged 
intensity maps and represent the 2, 3, 4, 5, 7 and 
9(7 levels with 1(7=4.4 mJy /beam. This rms level has 
been estimated from the thermal noise measured in 2005 
(22.8mJy/beam per channel; Table 1) divided by VN, 
where N is the number of channels. For the velocity 
range from 20 to 50kms~^, the number of channels is 
~27. As mentioned in Section 3, we use the 2005 ther- 
mal noise as a reference, since this is the largest r.m.s. 
noise measured between both epochs (Table 1). 

Figure 1 shows that the SiO line emission at intermedi- 
ate velocities has experienced signiflcant changes in mor- 
phology and intensity since 2005. The SiO emission in 
the 2005 image is clearly restricted to one condensation 
close to the L1448-mm protostar at offset (2,-2), emis- 
sion peak A, which is coincident with the molecular bul- 
let Rl (Guilloteau et al. 1992; Dutrey et al. 1997; Girart 
& Acord 2001). This emission is also seen in the higher- 
excitation transition SiO J~8^>-7 imaged by Hirano et al. 
(2010). In the 2009 map, however, SiO J=2^1 is also 
detected toward three other fainter clumps (seen at the 
2-4(7 level) along the direction of the high- velocity jet. 
These clumps are located at offsets (3,-9) for B, (6,-16) 
for C, and (9,-21) for D, and do not seem to have any 
counterpart in previous observations of the SiO J=2— >1 
emission carried out with the IRAM PdBI with compa- 
rable sensitivity (see Figures 2, 3 and 4 in Dutrey et al. 
1997). Indeed, the position- velocity diagrams of SiO ob- 
tained along the outflow (Dutrey et al. 1997) show that 
the region where the new B, C and D clumps are found 
(from 10" to 30" south the L1448-mm central source) was 
devoid of SiO emission in 1993, supporting the idea of 
an increase in the SiO intensity since 2005. Outside the 
primary beam, another very faint condensation, clump 
E at offset (11, -34), is also found along the jet prob- 
ably belonging to the same global structure of the new 
SiO clumps. This clump could be associated with the 
molecular bullet R2 already detected in the SiO maps of 
Dutrey et al. (1997). We note that the condensations B, 
C, D and E are unresolved in the direction perpendicular 
to the high- velocity jet with measured sizes of <4" . The 
derived deconvolved sizes are of <1.8". From all this, 
we propose that the new SiO condensations detected at 
intermediate velocities are due to an increase in the SiO 
gas-phase abundance generated by the recent interaction 
of young C-shocks with the ambient gas. 

To show that the changes in the morphology and in- 
tensity of the SiO emission found in the 2009 images are 
not a consequence of the non-identical UV coverage of the 
2005 and 2009 observations, in the left and middle panels 
of Figure 2 we present the dirty beams of these observa- 
tions for the region reported in Figure 1. In the right 
panel of Figure 2, we also show the subtraction of the 
two dirty beams, 9 pw h m C^OO^yO pw h m C^OOQ) ■ Dashed 
boxes show the regions where the new SiO condensations 
have been detected (see Figure 1). From Figure 2, it is 
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Table 1. Median values of the visibility noise per channel for all 
derived uv tables 





Chan. Res. 


Epoch 


VisibiUty 


Noise 




Bandwidth 


I'hcrmal 


Heal 


Imaginary 




80 MHz 
80 MHz 


312 kHz 
312 kHz 


2005 
2009 


22.8 
12.6 


23.1 
14.4 


23.2 

14.4 


mjy/beam 
mjy/beam 


20 MHz 
20 MHz 


39 kHz 
39 kHz 


2005 
2009 


21.8 

18.7 


22.0 
21.1 


21.7 
21.1 


mjy/beam 
mjy/beam 



clear that both beams are very similar with intensity dif- 
ferences of less than 4.5 mJy/beam (i.e. 0.045 after being 
normalized to the beam peak flux of 100 mJy/beani) for 
condensations A, B and C. This intensity difference is 
within the r.m.s. of our observations (4.4mJy/beani), 
which implies that the SiO intensity changes observed 
toward A, B, and C are likely due to time variability. 

Toward condensation D, the SiO emission reaches the 
4(7 detection level (i.e. 17.6mJy/beam). However, the 
absolute value of the intensity differences between the 
2005 and 2009 beams gets as large as 23mJy/beam, 
making the variable nature of the SiO emission toward 
this condensation uncertain. For clump E, the inten- 
sity differences between the beams reach 27mJy/beam 
in absolute value toward the north-west region of box 
E, where no variation in the SiO emission has been de- 
tected. In the south-east region of this condensation, 
the beam differences are less than 3 mJy/beam imply- 
ing that, although faint (intensity at the 2a level, i.e. 
8.8mJy/beam), the new SiO clump seen at the border 
of the primary beam could have also appeared recently. 

The increase of the SiO emission flux at intermediate 
velocities in 2009 is also evident from the right panels 
of Figure 1 , where we compare the SiO spectra averaged 
within condensations A, B, C, D and E for the 2005 and 
2009 observations, and smoothed to a velocity resolution 
of 2.2 km s~^. The intensity of the spatially averaged SiO 
spectra in Figure 1 are given in units of Jy. In Table 2 
(Columns 3 to 5), we report the SiO integrated inten- 
sity flux measured in both epochs toward A, B, C, D 
and E for the 20-50 km s~^ velocity range. In Table 2, we 
also include the SiO fluxes derived for the velocity ranges 
from -10 to 20kms^^ and from 50 to 80kms~^. The 
errors in the integrated intensity flux are calculated as 
(TArea—<^v5v^v, where (J is the rms noise level in the SiO 
spectra per channel (0.008 Jy and 0.006 Jy for the 2005 
and 2009 SiO spectra, respectively), 5v is the velocity 
resolution (2.2 kms^^), and Aw is the total width con- 
sidered in every velocity range (30kms~^). In Columns 
6 to 8 of Table 2, we show the ratios between the 2005 
and 2009 SiO fluxes, F(2009)/F(2005), and in Columns 
9 to 11, the absolute differences between these fluxes 
|F(2009)-F(2005)|. The uncertainties of these quantities 
have been estimated by propagating errors. 

From Figure 1 and Table 2, we find that the SiO line 
profiles have experienced significant flux variations since 
2005 in the 20-50 km s"-'^ velocity interval. Indeed, the 
measured SiO flux is factors 1.4-3 larger in 2009 for con- 
densations A, B, C, D and E. The absolute differences in 
the measured SiO flux range from 0.25 Jy kms^^ toward 
E to ~0.8 Jy kms^^ toward A, i.e. factors of ^2-7 larger 
than the estimated error of ±0.12 Jy kms^^. Although 



falling outside the primary beam of our observations, the 
averaged SiO line profiles toward E show an evolutionary 
trend for SiO to show slightly brighter emission at higher 
velocities in 2009 than in 2005. 

One may think that the increase in the SiO flux mea- 
sured in 2009 with respect to 2005 could be produced 
by an effect of differences in the data calibration process 
between both epochs. To check for this potential prob- 
lem, in Figure 3 we compare the images of the contin- 
uum emission at 87 GHz obtained in the 2005 and 2009 
observations. If calibration effects were responsible for 
the SiO flux increase, they should be apparent in the 
continuum images since the continuum emission flux is 
expected to remain constant with time. From Figure 3, 
we find that the peak continuum flux for the 2009 epoch 
is weaker (~9 mJy/beam) than in the 2005 observations 
( 15mJy/beam). These peak continuum fluxes are con- 
sistent with those previously measured by Guilloteau et 
al. (1992) also with the PdBI. Like in the SiO line im- 
ages, the differences between the 2005 and 2009 beams 
(right panel in Figure 3) remain below (<1.4mJy/beam) 
the rms noise level in the images (2 mJy/beam in 2005 
and 1.5mJy/beam in 2009). From all this, it is clear 
that the trend for the SiO line flux to increase from 2005 
to 2009 is the opposite to that found for the continuum 
emission between both epochs, and therefore cannot be 
attributed to calibration effects. We note that the dif- 
ference between the continuum peak fiuxes measured in 
2005 and 2009 (^6 mJy/beam) is significantly above the 
rms noise level of the 2005 image, but still within 3cr. 

From Table 2, we also find that the absolute differences 
for the measured SiO flux between 50 and 80kms~^ 
are similar to those for the velocity range from 30 to 
50kms^^. Indeed, small deviations in the SiO line pro- 
files can also be seen for the 50-80 km s~^ velocity interval 
(Figure 1). As shown in Figure 4, these deviations do not 
translate into significant variations in the morphology of 
the SiO emission toward clumps A and B. Toward C and 

D, slightly brighter SiO emission (at the 2-3(7 intensity 
level) is seen arising from these condensations. However, 
as mentioned before, the variable nature of this emission 
toward D is unclear, since the derived differences between 
the 2005 and 2009 beams are larger than 20mJy/bcam 
(i.e. than the 4(7 level in the 2005 and 2009 images). 

We note that the SiO fiux changes in the SiO emission 
for the intermediate- and high-velocity shock regimes to- 
ward condensations A, B and C, and probably toward 

E, are likely to be real because very small deviations are 
found for the low- velocity regime from -10 to 20kms^^. 
If these variations were produced by a low-level contribu- 
tion due to small differences in the UV coverage between 
the 2005 and 2009 observations, they would appear more 
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Fig. 1. — Left and middle panels: Averaged intensity images of the SiO J=2— ^-l emission from 20 to 50kms ^ observed toward L1448-mm 
in 2005 and 2009. Contour levels in both maps are 8.8 {2a), 13.2, 17.6, 22, 30.8 and 39.6 mjy/beam. The la level (4.4mjy/beam) is 
calculated from the thermal noise measured in the 2005 UV dataset (Section 4 and Tablet). Negative contours correspond to the —3a 
noise level (dotted contours). Arrows show the direction of the SiO high-velocity jet with velocities from 50 to 80kms~^. Filled star gives 
the location of the L1448-mm source (Maury et al. 2010). Solid ellipse shows the primary beam of our PdBI observations. Dashed boxes 
(A through E) delineate the regions where the 2005 and 2009 SiO spectra have been averaged. Beams are shown at the lower right corners 
and arc 3.80" x3. 35", P.A.=179° for the 2005 images, and 3.87" x3. 31", P.A.=157° for the 2009 maps. Right panels: SiO spectra averaged 
within condensations A, B, C, D and E for the 2005 (black histograms) and 2009 (red histograms) epochs. The units of the spatially 
averaged SiO spectra are given in units of Jy. 



dramatically at all velocities which clearly contrasts with 
our results shown in Table 2. 

4.2. Variability of the narrow SiO emission at the 
magnetic precursor components 

Figures 5 and 6 show the averaged intensity images of 
the SiO J— 2^1 emission observed in the two epochs 
for four different narrow SiO velocity components de- 
tected between 4.7kms~^ (the ambient cloud velocity) 
and lOkms"^. These components, Ml, M2, M3 and 
M4, have very narrow line profiles (linewidths of ^0.25- 
0.7kms~^), and are centered at radial velocities 5.2, 6.5, 
7.1 and 9.2kms~^, respectively. The rms noise level in 
every image has been calculated from the thermal noise 
shown in Table 1 divided by ^/N with N the number of 
channels used to generate the averaged intensity map. 
For the Ml and M3 components, we consider as refer- 
ence the rms noise measured in the 2005 observations 
(i.e. 7.7mJy/beam for Ml and 12.6mJy/beam for M3). 
For the M2 and M4 components, in general fainter than 
the Ml and M3 narrow emission, we have used the indi- 
vidual rms noise level obtained in every image (for M2, 
8.9mJy/beam and 7.6mJy/beam in the 2005 and 2009 



images, respectively; and for M4, 15.4 mJy /beam and 
13.2mJy/beam in the 2005 and 2009 maps). 

From Figure 5, we find that the narrow SiO emission 
at 5.2 km s"'^ (i.e. the Ml component) with a linewidth 
of only ^0.4 km s""'^, is the only component clearly de- 
tected in both epochs. This emission, previously re- 
ported by Jimenez-Serra et al. (2004) from single-dish 
IRAM 30 m observations, was interpreted as the signa- 
ture of the magnetic precursor of young C-shocks toward 
the L1448-mm outflow. The Ml emission peaks at (-2,- 
7), and although compact, shows an elongation in the 
direction perpendicular to the high-velocity jet (size of 
-4"x6", 0.005 pcx0.007pc). We note that this emission 
falls to the north of the secondary CO J=3— >-2 outflow 
reported by Hirano et al. (2010) and powered by the 
south-eastern source L1448-C(S) (Jorgensen et al. 2006). 
In addition, the narrow SiO kinematics do not exactly 
match those of the CO J=3— >-2 emission, suggesting that 
the Ml emission is not associated with this secondary 
outflow. 

From the averaged SiO spectra of Figure 5 (right pan- 
els), the derived SiO integrated flux for Ml in the 2009 
spectrum is slightly brighter than in the 2005 data by 
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Fig. 2. — Comparison of the beams obtained in tlie 2005 and 2009 observations with the PdBI (left and middle panels), and subtraction 
of the 2009 beam pattern from the 2005 beam (right panel). First (solid) contour and step level correspond to 10% the peak value of the 
beam. Negative (dotted) contours are -10%, -20% and -30% this maximum value. Color scale is shown on the right panel, and beam sizes 
are indicated at the upper part of the left and middle panels. Solid ellipse shows the primary beam (FOV~58"). Dashed boxes (right 
panel) delineate the regions where variability of the SiO emission has been detected (see Figure 1). 
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Fig. 3. — Images of the continuum emission at 87 GHz observed with the PdBI toward L1448-mm in the 2005 (left panel) and 2009 
(middle panel) epochs. The first contour and step levels are 6 (Scr) and 2mjy/beani in the 2005 image, and 4.5 (3(t) and 1.5nijy/beam 
in the 2009 map. The peak intensity of the continuum emission measured in 2005 was '~15mjy/beam, while in 2009 this intensity was 
~9mjy/beam. The synthesized beams are shown in the lower left corner of each panel and are 3.48" X 3.33", P.A.= 28° for the 2005 epoch, 
and 3.44" x3.37", P. A. =129° for the 2009 map. In the right panel of this Figure, we also show the subtraction of the 2009 beam from the 
2005 beam. Color scale and contour levels are as in Figure 2. 



20%±13%. The error in the relative difference of the 
SiO flux is large. However, as shown in Section 4.1, the 
flux increase by 20% should be considered as a lower limit 
since the continuum peak intensity in the 2009 images is 
weaker than in 2005. Therefore, it is possible that the 
narrow SiO emission arising from Ml is also variable. 

In contrast with Ml, the other three narrow SiO com- 
ponents (M2, M3 and M4) seem to be associated with 



transient structures that appear either in the 2005 or 
2009 images. For instance, for the M2 narrow emission 
at 6.5kms~^, one condensation is found toward clump 
B in the 2005 map, while two new M2 condensations ap- 
pear associated with clumps C and with the south-east 
edge of condensation D. From the SiO spectra averaged 
over these new condensations (Figure 5), we find that the 
SiO line profiles seem to have evolved since 2005 because 
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Fig. 4. — As in Figure 1, but for the velocity range from 50 to SOkms^"*^ (see vertical dotted lines in the right panels of this Figure). 
Contour levels in both maps are 13.2 (Scr), 17.6, 22, 30.8 and 39.6mjy/beam, with lcr=4.4 mjy/beam (Section 4 and Table 1). Negative 
contours correspond to the —3(7 noise level (dotted contours). 



their emission peaks appear more red-shifted within the 
shock, and their line wings are slightly brighter at veloc- 
ities ranging from 7 to lOkms"^. In addition, the M2 
condensation toward D is particularly interesting because 
it also shows an elongation in the direction perpendicular 
to the jet, as if tracing the interface between the recently 
shocked gas of clumps D and E. 

We note that the detection of the M2 emission com- 
ponents in 2009 could be a direct consequence of the 
higher sensitivity of the 2009 images. However, no simi- 
lar structure is seen in the 2005 map above the 2a level 
(17.8mJy/beam; left panel of Figure 5), which is below 
the 3(7 noise level (of 22.8 mJy) of the 2009 image. This 
suggests that the M2 narrow SiO lines, with linewidths 
of 0.7 km s^^, were likely not present in 2005. 

The M3 and M4 narrow SiO components (see Figure 6) 
have linewidths of ^0.25-0.6 km s~^ and their morphol- 
ogy is more compact (sizes of <4"). The MS and M4 
clumps are found toward A and B, likely forming part 
of a larger C-shock structure (left panels of Figure 6). 
The M3 emission seen at the interface between clumps 
D and E in 2005 seems to have vanished in the 2009 im- 
age, revealing the transient nature of the precursor phe- 
nomenon. It is remarkable that the M4 clump close to the 
L1448-mm protostar is relatively bright (peak intensity 
of ~0.05 Jy/beam) but remains unresolved in the 2009 
images. This implies that the magnetic precursor region 



toward this component is smaller than ~3.5" (0.004 pc). 

5. ORIGIN OF THE VARIABILITY OF THE THERMAL SIO 
LINE EMISSION TOWARD L1448-MM 

Variability of the atomic/molecular emission toward 
molecular outflows have usually focused on the detec- 
tion of proper motions associated with the high- velocity 
shocked gas in these objects (Reipurth & Bally 2001; 
Girart & Acord 2001; McGroarty, Ray & Froebrich 
2007). However, it has usually been neglected that the 
moderate- and intermediate-velocity shocked gas gener- 
ated in C-shocks, could also be subject to changes in the 
line emission arising from these regions. 

Our data show that the SiO emission at intermediate 
velocities toward the red-shifted lobe of the L1448-mm 
outflow has changed in morphology, line intensity, and 
line shape in only 4yrs. This is particularly true for con- 
densations A, B and C, where their measured SiO flux 
increase by factors 1.4-3 cannot be attributed to system- 
atic differences between the beams of the observations 
(Section4.1). This flux increase suggests that new ma- 
terial has likely entered the shock, enhancing the gas- 
phase abundance of SiO toward these regions. The lack 
of detection of SiO in previous interferometric observa- 
tions toward the location of the new SiO condensations 
(Dutrcy et al. 1997), also favours the idea of a recent SiO 
enhancement toward the red-shifted lobe of the L1448- 
mm outflow. Although not predicted by the models of 
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Fig. 5. — Averaged intensity maps of the SiO J=2— >1 emission observed toward L1448-mm in 2005 (left panel) and 2009 (central 
panel) for the magnetic precursor components at 5.2kms~^ (Ml; red contours) and at 6.5kms~^ (M2; blue contours). These images are 
superimposed on the maps of the intermediate- velocity SiO emission shown in Figure 1 (grey scale and black contours). The first contour 
and step levels for the Ml component are 23.1 (Scr) and 23.1mjy/beam in the 2005 and 2009 maps. The Itr level (7.7mjy/beam) has 
been calculated from the thermal noise of the 2005 data (Table 1). For the M2 component, the first contour and step levels for the 2005 
image are 17.8 (2(t) and 8.9mjy/beam, and for the 2009 map, 15.2 (2ct) and 7.6mjy/beam. The individual rms for every map has been 
estimated from the thermal noise of Tablet. Filled star and arrows are as in Figure 1. Beams are shown at the lower left corners. Thick 
ellipse delineate the primary beam of the PdBI observations. Dashed boxes show the regions where the SiO spectra have been averaged 
(left panels). Dotted vertical lines in the spectra indicate the central radial velocities of the Ml and M2 precursor components. 



Table 2. Ratios and absolute differences in the measured SiO flux 
between the 2005 and 2009 spectra. 



F [Jykms-1] F(2009)/F(2005) 

Epoch -10^20 20^50 50^80 -10^20 20^50 50^80 

(kms~^) (kms~^) (kms~^) (kms~^) (kms"-'-) (kms"-"-) 



|F(2009)-F(2005)| [Jykms-i] 

-10^20 20^50 50^80 

(kms^-*-) (kms"-"-) (kms"-*-) 



D 



E 



2005 



<0.2=' 



1.93(7)'' 5.97(7) 



2009 


<0.15 


2.69(5) 


6.60(5) 


2005 


0.24(7) 


<0.2 


3.33(7) 


2009 


<0.15 


0.65(5) 


3.63(5) 


2005 


<0.2 


0.23(7) 


0.63(7) 


2009 


<0.15 


0.71(5) 


1.00(5) 


2005 


<0.2 


0.25(7) 


0.31(7) 


2009 


0.16(5) 


0.68(5) 


0.80(5) 


2005 


<0.2 


<0.2 


<0.2 


2009 


<0.15 


0.44(5) 


0.37(5) 



1.40±0.05 1.106±0.015 

<0.63 >3.3±1.1 1.09±0.03 <0.09 

3.1±0.9 1.6±0.2 

>0.8 2.7±0.7 2.6±0.6 <0.04 

>2.3±0.8 >1.9±0.7 



0.76±0.12 0.63±0.12 

0.45±0.12 0.29±0.12 

0.47±0.12 0.37±0.12 

0.43±0.12 0.49±0.12 

0.25±0.12 0.18±0.12 



^Upper limits correspond to the SaAma level, where a Area is calculated as cry/SvAv (see Section 4.1). 
''Uncertainties in the SiO integrated intensity flux correspond to the (TArea level. 



Jimenez-Serra et al. (2009) for such short time-scales, 
the SiO line profiles could also be variable at high veloci- 
ties with a moderate increase in their integrated intensity 
flux. 

The origin for such a coherent variability in the SiO 
emission over spatial scales of 750-1000 AU (i.e. over 
the 30-40" mapped with the PdBI toward the red-shifted 
lobe of the L1448-mm outflow) could be related to the 
formation of C-shocks in the internal working surface 
between the high-velocity jet and the ambient mate- 



rial. These shocks would tend to propagate sideways 
and would participate in the expansion processes associ- 
ated with the high- velocity jet (see e.g. Cabrit 1995). As 
shown by the extended SiO and H2 shock-excited emis- 
sion reported by Dutrey et al. (1997) and Neufeld et al. 
(2009), the L1448-mm jet has already travelled along dis- 
tances larger than >2', explaining the coverage of at least 
30-40" where variability of the SiO emission has been 
detected. In Section 4.1, we have shown that the decon- 
volved size of the B, C, D and E condensations in the 
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Fig. 6. — Same as in Figured, but for the precursor components M3, at 7.1kms~^ (red contours), and M4, at 9.2kms~^ (blue contours). 
The first contour and step levels for the M3 component are 25.2 {2a-) and 12.6 mjy/beam in the 2005 and 2009 maps. For the M4 component, 
the first contour and step levels are 30.8 (2cr) and 15.4mjy/beam in the 2005 image, and 26.4 {2cr) and 13.2mjy/beam for the 2009 map. 
Dashed boxes limit the region where we have averaged the SiO emission toward the M3 condensation. For the M4 component, we only 
show the SiO spectra extracted from the emission peaks of the M4 condensations (right panels). Dotted vertical lines in the SiO spectra 
indicate the radial central velocities for M3 and M4. 



direction perpendicular to the high- velocity jet is <1.8" 
(i.e. 0.002 pc at a distance of 250 pc). Assuming that 
C-shocks generated at the working surface travel sym- 
metricaly and sideways half that distance from the jet, 
we estimate a propagation velocity for those C-shocks of 
<240kms~^ in the 4yrs time-interval between the 2005 
and 2009 observations. These velocities are consistent 
with those derived by Girart & Acord (2001, of few hun- 
dred kms"^) from the proper motions of the SiO high- 
velocity bullets observed toward the red-shifted lobe of 
the L1448-mm outflow. 

We would like to stress that the presence of C-shocks 
with such high velocities could be possible within an sce- 
nario of a highly magnetized medium with a large pre- 
shock magnetic induction, B, transverse to the direction 
of the flow. Indeed, Le Bourlot et al. (2002) showed 
that, for large values of B (of the order of some mG), 
the critical velocities of C-shocks can be as high as few 
100 km s^^, which are similar to those derived from our 
observations. 

The variability of the thermal SiO emission toward 
L1448-mm is also supported by our high-spectral reso- 
lution SiO data, which show regions (the M2 precursor 
component) where the narrow SiO line profiles seem to 
have evolved since 2005. The emission peak of these lines 
have been red-shifted within the shock, and have devel- 
oped slightly brighter broad SiO emission as predicted 
by shock models (Jimenez-Serra et al. 2009). Although 
the broad SiO condensation E falls outside the primary 
beam of our observations, this trend in the evolution of 



the SiO line profile since 2005 is also found toward this 
condensation. 

The narrow SiO emission toward the Ml component 
also seems to have experienced a small increase in its 
integrated intensity flux (by 20±13%). This component, 
previously reported by single-dish observations (Jimenez- 
Serra et al. 2004), is known to be correlated with an 
enhancement and over-excitation of the ion fluid (rep- 
resented by the H^'^CO+ molecule) with respect to the 
neutral one (with molecules such as HN^'^C and H^'^CN; 
see e.g. Jimenez-Serra et al. 2006; Roberts et al. 2010), as 
expected within the magnetic precursor scenario. How- 
ever, although predicted by our C-shock modelling re- 
sults (Jimenez-Serra et al. 2009), no significant variation 
in the linewidth of narrow SiO is detected for the Ml 
component toward the L1448-mm outflow. 

The other narrow SiO components (M2, M3 and 
M4) are fainter and their linewidths are as narrow as 
~0.25kms~^. Even in the case that no turbulence was 
present in the medium, these linewidths would imply 
a kinetic temperature of <60K for the SiO gas, which 
could only be attained at the early precursor stage or 
at the evolved postshock gas (see e.g. Flower & Pineau 
des Forets 2003; Gusdorf et al. 2008). Since these com- 
ponents seem to be transient (Section 4.2), it is likely 
that they arise from the precursor region of the shock. 
Reflected shock waves within the postshock gas could 
not be responsible for the narrow SiO lines in this out- 
flow, because their ion-neutral velocity drift never ex- 
ceeds 6kms~^ (Ashmore et al. 2010). This is the lower 
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limit drift velocity required to sputter detectable abun- 
dances of SiO in the gas-phase, from the mantles of dust 
grains (Caselli et al. 1997; Jimenez-Serra et al. 2008, 
2009). 

The morphology of the M2 emitting regions appears 
elongated in the direction perpendicular to the high- 
velocity jet, resembling the bow-shock structures ob- 
served toward other jets and outflows (e.g. as in HHl-2 
or Orion; Eisloffcl, Mundt, & Bohm 1994; Davis et al. 
2009). In addition, these regions remain practically un- 
resolved in the direction of the jet (<4", 1.5-1.8x10^® cm; 
Section 4.2), and are generally found at the interface be- 
tween the recently shocked condensations A, B, C, D 
and E detected at intermediate velocities in the 2009 
maps. This indicates that the interface/precursor regions 
typically have very small length scales as predicted by 
C-shock modeling (~3-20xl0i5cm; Gusdorf et al. 2008; 
Jimenez-Serra et al. 2008). This is also consistent with 
the fact that the M4 narrow SiO line arises from unre- 
solved regions (Section4.2). 

We note that, although we compare our results with 
those from steady-state C-shock models, recent modeling 
of the transient evolution of C-shocks predict the pres- 
ence of a magnetic precursor within the physical struc- 
ture of the shock, as a consequence of its propagation 
through a decreasing density medium such as that found 
in molecular outflows (Ashmore et al. 2010). 

6. CONCLUSIONS 

In this paper, we present interferometric images of the 
thermal emission of the SiO J=2^1 rotational line ob- 
served toward the L1448-mm outflow at two different 
epochs (2005 and 2009). The PdBI images reveal that 
this emission is variable in only 4yrs. Not only the mor- 
phology of the SiO emission, but also its line profiles. 



have significantly changed since 2005, with a clear in- 
crease of the SiO flux at intermediate- velocities, and pos- 
sibly at high-velocities. This suggests that new material 
has recently entered the shock. 

For the precursor component, we have detected sev- 
eral narrow SiO velocity components (Ml, M2, M3 and 
M4) with an elongated morphology in the direction per- 
pendicular to the high-velocity jet. These features re- 
semble those found in bow-shocks. Although these com- 
ponents could have appeared recently, new multi-epoch 
and higher sensitivity interferometric images are needed 
to clearly establish the variable nature of the narrow SiO 
emission detected toward L1448-mm. 

We speculate that the variability of the thermal SiO 
J=2^1 emission observed toward the L1448-mm outflow 
could be explained by the time-dependent propagation of 
very young C-shocks (with their precursors) through the 
ambient molecular cloud. Monitoring of the variability 
of the thermal SiO emission toward young molecular out- 
flows could be a useful tool to study the evolution of these 
objects, and their impact on the surrounding medium. 
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